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Uncertainty  of  the  NIST 
Electrooptic  Sampling  System 


Dylan  Williams 
Electromagneucs  Division 

Paul  Hale  and  Tracy  Clement 
Optoelectronics  Division 

Chih-^Eng  Wang 
Statistical  Engineering  Division 

National  Institute  of  Standairds  and  Technology 

We  analyze  the  uncertainU'  of  XLST's  electrooptic  sampling  system.  The  system  measures  voltage 
waveforms  m  a  coplanar  waveguide  fabricated  on  an  electrooptic  LiTaO-  wafer.  We  use  the  system 
to  measure  the  voltage  waveform  injected  by  a  photodetector  through  a  microwave  probe  into  this 
coplanar  waveguide.  We  then  determine  the  voltage  the  photodetector  would  suppl\-  to  a  50  Q  load 
at  its  coaxial  connector.  We  calculate  this  voltage  using  an  electrical  mismatch  correction  that 
accounts  for  the  effects  of  the  probe  and  coplanar  wa\'eguide  on  the  measured  waveform. 

Ke^'vvords:  electrical  phase:  electrooptic  sampling;  mismatch  correction:  photodetector  calibration: 
photodiode  calibration. 

1.  Introduction 

We  present  an  unceitaint}'  anah  sis  for  measurements  performed  with  the  NTST  electrooptic  sampling 
(EOS)  system  described  in  Refs.  [1]  and  [2].  The  EOS  system  measures  the  voltage  waveform 
injected  by  a  photodetector  into  a  coplanar  wa\eguide  i  CPWi  fabricated  on  an  electrooptic  LiTaO- 
wafer.  During  post  processing,  we  take  a  Fourier  transform  of  the  temporal  waveform  measured  in 
the  CPW.  Using  the  mismatch  corrections  described  in  Ref.  [IT  we  then  determine  the  magnitude 
and  phase  of  the  frequency  response  of  the  \oltage  the  photodetector  would  deli\er  to  a  50  Q  load 
at  its  coaxial  porL 

We  currentiy  perform  r.vo  upes  of  micasurements  on  the  EOS  s>-stem:  measurements  of 
photodetectors  prodded  by  customers  and  measurements  of  our  own  check  standards.  When  we 
measure  a  customer's  photodetector.  we  repon  the  mean  of  n.  =  3  measurements.  \Mien  we  report 
a  measurement  of  our  check  standard,  we  average  sets  of  three  measurements.  In  some  cases, 
these  additional  measurements  will  reduce  our  measurement  uncertainties  significanth'.  Some 
customers  use  the  NIST  check  standard  to  cahbrate  their  instruments  directiy. 

In  the  followine.  we  separate!}"  estimate  systematic  p.pe-B  and  statistically  derived  t>pe-A 
uncertainties  [3].  We  also  separately  estimate  the  components  of  our  r.pe-A  uncertainr."  due  to 


1 


repeatability  within  each  set  of  3  measurements,  and  the  long-term  reproducibility  between  sets 
of  three  measurements.  We  also  use  a  magnitude/phase  representation  for  the  data,  rather  than  a 
real/imaginary  representation.  This  is  reasonable  because  the  magnitude  is  much  larger  than  the 
square  root  of  its  variance  s^,  and  the  variance  Sq  of  the  phase  is  small.  Also,  throughout  the  analysis 
we  treat  the  magnitude  and  phase  separately. 

1.1  Check-standard  Uncertainty 

We  report  the  mean  of  sets  of  n^=3  measurements  when  characterizing  the  check  standard.  To 
estimate  the  component  of  our  type-A  uncertainty  in  these  measurements  due  to  short-term 
repeatability,  we  average  the  variances  from  each  set  of  measurements.  We  also  separately 
estimate  the  component  due  to  long-term  reproducibility  using  a  variance-component  model.  Finally, 
we  use  a  Monte-Carlo  simulator  to  estimate  the  uncertainty  in  the  measurement  due  to  systematic 
errors,  which  in  our  case  are  all  of  type  B. 

1.2  Uncertainty  in  Measurements  of  a  Customer-supplied  Photodetector 

To  estimate  the  component  of  the  type-A  uncertainty  in  the  measurements  of  a  customer's 
photodetector  due  to  short-term  repeatability,  we  perform  each  measurement  of  the  photodetector' s 
frequency  response  three  times  (n^  =  3).  We  perform  only  three  measurements  because  each 
measurement  takes  approximately  100  hours  to  perform,  and  time  constraints  limit  us  to  three  repeat 
measurements  per  photodetector. 

We  report  the  mean  of  these  three  measurements  of  the  magnitude  and  phase  of  the  photodetector' s 
frequency  response  when  characterizing  a  customer's  photodetector.  We  rely  on  our  check-standard 
measurements  to  estimate  the  component  of  our  type-A  uncertainty  in  our  measurements  due  to  long- 
term  reproducibility.  Finally,  we  use  a  Monte-Carlo  simulator  to  estimate  the  uncertainty  in  the 
measurement  due  to  systematic  errors. 

1.3  Overview  of  the  Uncertainty  Analysis 

For  both  our  customer  and  check-standard  measurements,  we  combine  the  repeatability, 
reproducibility,  and  systematic  components  of  uncertainty  to  estimate  the  combined  uncertainty  and 
confidence  intervals  for  our  measurements.  Throughout,  we  add  a  subscript  'r'  to  quantities 
associated  with  components  of  uncertainty  due  to  imperfect  short-term  repeatability  of  measurements 
within  a  single  set,  a  subscript  'R'  to  quantities  associated  with  the  component  of  our  uncertainty  due 
to  imperfect  long-term  reproducibiUty  between  sets,  and  a  subscript 's'  to  quantities  associated  with 
our  type-B  uncertainties  caused  by  the  systematic  errors  we  identified.  Due  to  the  small  number  of 
measurements  we  perform,  we  do  not  estimate  the  correlation  between  the  magnitude  and  phase 
measurements. 

We  combine  these  uncertainties  following  the  recommendations  of  Ref.  [3]  in  four  consecutive 
steps.  These  are: 
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(1)  Calculate  the  sample  variances  of  the  systematic  (type-B),  repeatability  (type-A),  and  long- 
term  reproducibility  (type-A)  errors  separately.' 

(2)  Estimate  the  standard  uncertainties  u  in  either  the  quantity  of  interest  or  its  mean,  as  appropriate, 
from  these  variances.'^ 

(3)  Calculate  the  combined  uncertainty  from  the  sum  of  squares  of  the  standard  uncertainties 
associated  with  the  systematic,  repeatability,  and  reproducibility  errors.  The  square  of  the 
combined  uncertainty  is  an  estimate  of  the  variance  of  the  averaged  magnitude  and  phase 
responses  we  report. 

(4)  Estimate  the  coverage  idLCiorh^^  and  expanded  uncertainty  U  corresponding  to  95  %  confidence 
intervals  for  the  means  we  report  from  the  combined  uncertainty  and  the  number  of  degrees  of 
freedom  v  associated  with  each  of  the  standard  uncertainties  u.  The  number  of  degrees  of 
freedom  v  reflects  the  confidence  we  have  in  the  standard  uncertainties  we  use  to  form  the 
combined  uncertainty.  The  coverage  factor  ^5  is  the  factor  by  which  we  must  multiply  the 
combined  uncertainty  to  form  the  expanded  uncertainty  U  and  95  %  confidence  intervals.  The 
coverage  factor  k^^  depends  on  the  number  of  effective  degrees  of  freedom  Vgf^,  which  in  turn 
depends  on  both  the  relative  sizes  of  each  component  of  the  combined  uncertainty  and  the 
number  of  degrees  of  freedom  of  each  component. 

2.     Uncertainty  due  to  Systematic  Sources  of  Error 

We  estimate  our  systematic  errors  in  the  three  magnitude-response  measurements  m,  and  the  three 
phase-response  measurements  6,  of  the  customer's  photodetector  with  a  Monte-Carlo  simulator.  For 
each  of  the  three  measurements,  we  run  our  Monte-Carlo  simulator  first  with  no  errors,  then  =  100 
additional  times  with  the  systematic  errors  detailed  in  the  rest  of  this  report  added  into  the 
measurements.  At  each  frequency  point  we  calculate  the  100  magnitude  and  phase  responses  of  the 
photodetector  in  the  presence  of  the  simulated  errors.  Finally,  using  =  100,  we  determine  the 
sample  variances  s^  jj'  and  s^  q^  of  the  100  simulated  values  of  the  magnitude  and  phase  for  each  of 
the  three  measurements  wi,-  and  0,.'  These  variances  estimate  the  standard  uncertainty  in  the  m,  and 
0,  due  to  systematic  sources  that  do  not  change  from  measurement  to  measurement,  and  are  all  of 
typeB. 

Since  the  sources  of  the  systematic  errors  are  similar  in  all  of  our  measurements,  and  the 
measurements  themselves  are  quite  close,  we  have  5^  ,,,1^  =  s^^^2  ~  ^s.m^  ^"^^ f  s.ei"^  ~  ^^m  ~  Thus 
we  estimate  our  standard  uncertainties  u^^  and     g  in  the  means  m  and  0  of  the  m,  and  0,  due  to 


'We  calculate  sample  variances  s  from  s'^  =  iyryf^  where  n  is  the  number  of  samples,  the    are  the 

individual  samples,  and  y  is  the  mean  of  the  j,. 

^  We  estimate  the  standard  uncertainty  of  a  quantity  as  s,  the  square  root  of  its  variance.  We  estimate  uncertainty  of 

means  as  \js^ln,  where  n  is  the  number  of  samples.  (Note  that,  since  averaging  more  values  increases  confidence  in 
a  mean,  the  standard  uncertainty  of  the  mean  of  a  quantity  is  smaller  than  s,  the  square  root  of  the  quantity's 
variance.). 
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these  systematic  sources  of  error  from  the  average  of  the  variances  s^  „^^  and  s^  q,^  using 


5^  ^s,mi  5^  ^s,ei 


2  j=l  ,        2  ,  =  1  (1) 

"s.rii  =    and   w  Q  =    , 

where  n^,  the  number  of  measurements,  is  equal  to  three.  We  assigned  an  infinite  number  of  degrees 
of  freedom     to     -  and    s,  as  recommended  in  [3]. 

s         s.tn  s,u ' 

3.     Variance  due  to  Short-term  Measurement  Repeatability 

The  analysis  of  the  random  part  of  our  uncertainty  begins  with  an  analysis  of  the  component  due  to 
short-term  measurement  repeatability.  To  estimate  our  repeatability  error  within  each  set  of  three 
measurements,  we  estimate  the  sample  variances  s,  J^  and  ^^.e^  of  the  three  measurements  m,-  and  6, 
using 


\m  =   ;         and  J  Q  =   . 

4.     Variance  due  to  Long-term  Measurement  Reproducibility 

hi  addition  to  the  short-term  measurement  repeatability  within  a  set  of  three  measurements  discussed 
in  the  last  section,  our  measurements  within  different  sets  may  also  differ  due  to  long-term 
reproducibility.  We  also  add  the  uncertainty  due  to  long-term  measurement  reproducibility  to  our 

2  2 

overall  uncertainty.  We  determine  the  variances  and  s^-^  characterizing  the  component  of  long- 
term  reproducibility  by  repeating  sets  of  three  check-standard  measurements.  We  do  this  using 

-2  -2 

« 2  2  ^T,m  1-2  2  -^r.e 

=  ^R,S  -  —   and   s^^  =  s^  ^  -  —  ,  (3) 

—  2  — 

which  separates  out  the  components  of  varience  s^^  and        due  to  short-term  measurement 

2  2 

repeatability  from  the  total  variances       and       of  the  means  of  our  sets  of  three  measurements 

 2   2 

performed  over  a  longer  time  span.  We  estimate  the  average  sample  variances  s^^  and  s^q  in  Eq. 

2  2 

(3)  due  to  short-term  repeatability  error  by  averaging  the  s^^  and  the  s^q  over  the     sets  of  12^=  3 

2  2 

measurements.  If  either  ij^-  or       in  Eq.  (3)  is  negative,  we  set  it  to  zero. 
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5.     Combining  Uncertainties 


We  estimate  the  combined  uncertainties  «-  and  Uq  of  m  and  0  from  the  long-term  reproducibihty 
components  Wj^-  and  u-^q  of  uncertainty,  the  short-term  repeatability  components  w^,-  and  u^q  of 
uncertainty,  and  the  systematic  components      -  and    g  of  uncertainty  using 

«l  =         +  ^r,m  +  «s'^   and   w|  =        +         +  w^-  .  (4) 

Thus  the  combined  uncertainties  w-  and  Wg  include  both  measurement  repeatability,  reproducibility, 

and  systematic  sources  of  uncertainty:  they  represent  our  best  estimate  of  the  standard  uncertainty 
of  our  results.  The  values  of    -  and    g  are  determined  by  (1).  The  values  of       ,      ,  Wj^q  ,  and  q 

will  be  described  below. 

5.1.  Customer  Measurements 

As  we  have  only  just  recently  finished  our  first  long-term  reproducibility  study,  for  our  customer 
measurements  we  have  to  date  neglected  Wj^-  and  u-^q  in  Eq.  (4).  When  we  are  ready  to  include 

and  u^Q  in  Eq.  (4),  we  will  determine  them  from  our  check-standard  measurements  with 

^Im  =  s^im   and   w^g  =  s^Q  .  (5) 

The  number  of  degrees  of  freedom  Vj^  associated  with  Wj^-  and  u-^q  is  np,  -  1.  For  customer 
measurements,  we  determine  the  standard  uncertainties       and  u^q  from 

2  2 

"r,m  =           and  M  g  =  —  ,  (6) 

where  n^,  the  number  of  measurements  used  to  form  the  mean,  is  equal  to  three.  The  number  of 
degrees  of  freedom    associated  with  u^^  and  u^q  is  n,-  1. 

5.2.  Check-standard  Measurements 

Since  we  average  measurements  to  characterize  our  check-standards,  we  do  not  use  (5)  and  (6) 
to  determine  the  uncertainties  in  Eq.  (4).  For  our  check-standard  measurements,  we  determine  the 
standard  uncertainties     -  andw^  s  due  to  long-term  measurement  reproducibility  with 
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«2  -  2 

2  "^R^m  2  "^ILG 

"r,iR  =    and  =           .  (7) 


Again,  the  number  of  degrees  of  freedom  Vr  associated  with  u^^-  andwj^g  is  /7r  -  1.  Since  for  our 

check  standard  we  report  a  mean  of  a  number  of  measurement  sets,  we  assign  our  short-term 
uncertainty  with 

-2  -2 

u  -  =    and    M  Q  =    ,  (8) 

"r  "r  '          «R  «r 

rather  than  with  (6).  Here  the  number  of  degrees  of  freedom  associated  with  w^.-  and  u^q  is 
Mr  K  -  !)• 

6.     Expanded  Uncertainty  and  95  %  Confidence  Intervals 

To  determine  our  expanded  uncertainties  U^^^  and  U^^q  corresponding  to  95  %  confidence 

intervals  for  in  and  0,  we  must  take  into  account  our  relative  confidence  in  the  standard 
uncertainties  u^-,  u^q  ,       and  u^q.  This  requires  determining  the  "coverage  factor"     [3]  from 

the  combined  uncertainties  u-  and  Ms . 

m  ti 

We  determine  kg^  from  v^jf,  the  effective  number  of  degrees  of  freedom  in  the  measurements.  We 
determine  v^f^  from  estimates  of  the  number  of  degrees  of  freedom  of  our  systematic  errors  and 
the  number  of  degrees  of  freedom  and  Vr  of  our  repeatability  and  reproducibility  with  the  Welch- 
Satterthwaite  formula 


and 


"im        \m  Km 

— —  +  — —  +  — — 

^R 


4 
"9 


^    +    ^    H-    ^  ^^^^ 


as  recommended  in  Appendix  B.3  of  Ref.  [3].  We  calculate  our  coverage  factors  kg^  from  = 
^95(^eff)'  where  tg^iv^ff)  is  the  two-sided  95'*'  percentile  of  Student's  r-distribution  with  v^^f  degrees  of 
freedom.  Finally,  we  determine  our  expanded  uncertainties  of  in  and  0  from  U^^  -  =  kg^^  m-  and 
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Figure  1 .  Representati%-e  measurement  of  a  check- 
standard  voltage  spectral  der^ity. 

^5,9  =  ^5,6^3  ■  ^^re  IS  a  9; 
that  the  true  values  of  m  and  6  lie  wuhm  the 
expanded  uncertainty'  inter\-als  m=  L\.  -  and 


Representative  Measurement  Result 


Figure  2.  Normalized  energy'  spectral  density' 
corresponding  to  the  voltage  of  Fig.  1. 
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From  each  measurement  on  the  electrooptic  I 
sampHng  system  we  form  both  an  absolute  and  I 
a  relative  voltage  magnitude  response  of  the  g;-:2. 
photodetector.  and  a  linearized  phase  resDonse. 

We  represent  the  absolute  voltase  response  of  Representative  uncertamt>-  in  our  check-standard 

,      ,        ,  ,  ,  _  ,     measurements  ot  spectral  densit%". 

the  photodetector  as  the  spectral  density"  of  the 

voltage  generated  by  the  detector  normahzed  to  the  charge  drawn  by  the  detector  when  it  is  excited 
by  a  narrow  optical  pulse.  This  spectral  density  has  units  of  V/C-Hz.  and  corresponds  to  the  voltage 
output  that  the  detector  would  generate  across  a  50  Q  load  connected  at  its  coaxial  output  pon  o\"er 
a  bandwidth  of  one  hertz  after  being  excited  by  an  optical  impulse  that  drew  one  coulomb  of  charge 
through  the  photodetector"  s  electrical  bias  port.  This  is  useful  because  we  can  better  measure  the  bias 
current  generated  by  the  photodetector  m  response  to  an  optical  impulse  than  determine  the  power 
in  each  optical  impulse.  Figure  1  shows  a  representati\"e  measurement  of  the  spectral  density  of  a 
detector's  \'oltage  response. 


We  repon  our  measurements  on  a  200  MHz  grid  to  110  GHz.  We  extract  these  measurements  from 
a  time  record  that  is  only  about  2  ns  long,  which  we  pad  to  5  ns  before  calculating  the  Fourier 
transforms.  Because  the  time  record  is  short,  measured  results  represent  average  values  o\  er  a 
bandwidth  of  rough!}'  500  MHz.  Furthermore,  small  v  oltage  offsets  in  the  temporal  measurements 
manifest  themselves  as  energy  between  DC  and  500  MHz.  We  have  not  yet  fully  characterized  or 
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Figure  4.  Comparison  of  five  check-standard  phase 
measurements. 
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Figure.  5.  Representative  uncertainty  in  our  check-standard 
measurements  of  phase. 


developed  methods  to  eliminate  these  small  offsets,  so  our  measurements  and  error  estimates  may 
not  reliably  reflect  the  actual  distribution  of  power  below  about  500  MHz. 

We  also  measure  a  normalized  version  of  the  voltage  spectral  density.  The  quantity  we  determine 
from  each  measurement  is  the  energy  that  the  photodetector  generates  in  a  50  Q  load  in  a  1  hertz 
bandwidth  normalized  to  the  total  energy  we  measure  from  the  detector  over  our  total  measurement 
bandwidth  (usually  110  GHz).  Figure  2  shows  the  normalized  energy  density  corresponding  to  the 
measurement  in  Fig.  1  expressed  in  decibels.  The  values  plotted  in  the  figure  are  on  the  order  of  10" ' ' 
because  the  energy  in  a  1  Hz  bandwidth  is  small  compared  to  the  total  energy  in  the  110  GHz 
spectrum  we  measure.  The  relative  scale  on  the  right  facilitates  determining  quantities  such  as  the 
detector's  3  dB  bandwidth;  Fig.  3  shows  the  associated  uncertainty  in  the  measurement. 

We  also  measure  the  relative  phase  response  of  the  detectors.  We  report  the  phase  of  the  voltage  the 
detector  will  generate  across  a  50  Q  load  after  subtracting  a  time  delay  from  the  phase  response  that 
minimizes  the  measured  phase  below  30  GHz  in  the  least-squares  sense.  Figure  4  shows  the 
measured  phase  response  of  one  of  our  check  standards.  It  compares  5  sets  of  three  measurements 
to  the  mean  of  these  five  sets  of  measurements. 


We  report  the  standard  and  expanded  uncertainty  for  each  measurement  as  well.  Figure  5  plots  the 
uncertainty  in  phase  for  the  check-standard  measurements  presented  in  Fig.  4.  For  this  measurement, 
which  is  an  average  of  the  mean  of  five  measurement  sets,  the  systematic  errors  dominate.  Keep  in 
mind  that  these  results  are  only  representative.  Not  only  will  results  vary  from  photodetector  to 
photodetector,  but  the  uncertainties  we  achieve  depend  on  experimental  conditions  and  may  change 
from  measurement  to  measurement.  Furthermore,  we  are  not  able  to  characterize  the  response  of  a 
customer's  photodetector  with  the  same  accuracy  as  our  check  standard,  because  we  are  unable  to 
perform  more  than  three  measurements  on  each  customer's  photodetector. 
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Figure  6.  The  electrical  model  used  to  represent  the  errors  in  the  electroopdc  sampling  system. 


8.     Specific  Forms  Used  to  Represent  Systematic  Errors 


In  what  follows,  we  derive  the  form  of  each  component  of  systematic  error  we  consider  to  determine 
the  variances  s.^_^[  ands,  gf .  These  error  mechanisms  can  be  divided  into  two  categories,  those  having 
to  do  with  the  mismatch  correction  (SOLT  calibration  errors,  reflection  coefficient  measurement 
errors,  and  on-wafer  scattering  parameter  measurement  errors)  and  those  having  to  do  with  the  finite 
impulse  response  of  the  EOS  system  itself.  We  model  most  of  the  errors  having  to  do  with  the 
mismatch  corrections  as  two-port  electrical  circuit  elements  m  the  electrical  circuits  between  the 
photodetector  and  the  on-wafer  reference  plane.  We  model  the  errors  having  to  do  with  the  finite 
impulse  response  of  the  EOS  system  as  multiphcative  errors  in  the  frequency  domain. 


Figure  6  shows  the  electncal  model  we  use  to  represent  the  errors  due  to  imperfect  mismatch 
corrections  in  the  electrooptic  sampling  system.  The  voltage  represents  the  voltage  across  the 
detector,  and  the  voltage  V^^w  represents  the  voltage  that  the  electrooptic  sampling  system  measures 
in  the  c opl an ar  waveguide  (CPW)  transmission  line.  Parameters  denoted  by  a  F  represent  reflection 
coefficients,  those  by  an  5  represent  scattering-parameter  matrices,  and  those  by  a  T  represent 
cascade  matrices.  The  function  T(5)  in  the  equation  embedded  in  the  figure  represents  the  cascade 
matrix  corresponding  to  its  argument,  which  is  a  scattering-parameter  matrix  5.  F^'  represents  the 
reflection  coefficient  of  the  detector  and  its  measurement  errors  (see  section  8.4.2).  5solt  ^^e  errors 
in  the  coaxial  SOLT  calibration,  the  errors  in  the  on-wafer  TRL  calibration  m  the  coplanar 
waveguide  transmission  lines,  and  ^he  reflection  coefficient  of  the  CPW  resistor.  The  quantities 
Tp,  5i,  ,  SJ\  and  AS'  represent  errors  in  the  on-wafer  TRL  calibration,  and  are  discussed  in  detail 
later. 


The  equations  in  the  rest  of  this  document  descnbe  m  greater  detail  how  we  add  systematic  errors 
into  the  Monte-Carlo  simulations  we  use  to  determine  the  variances  s^_^^'  and  s^Qf.  We  use  the 
s\mbol  \i  to  denote  random  variables  with  a  zero  mean  and  Gaussian  distribution  having  a  variance 
of  L  Using  this  notation,  we  would  use  m  +  o\ilo  represent  a  random  vanable  having  a  Gaussian 
distribution  with  mean  m  and  variance  o-  (i.e.,  a  Gaussian  distnbution  having  mean  m  and  standard 
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deviation  a). 


The  symbol  \x  with  no  argument  denotes  a  single  variable  chosen  once  randomly  during  each 
iteration  of  the  Monte-Carlo  simulator.  The  symbol  denotes  a  variable  chosen  randomly  during 
each  iteration  at  each  frequency  /.  No  |i  in  any  equation  is  correlated  to  any  other  |i  in  any  other 
equation.  When  we  use  more  than  one  random  variable  within  the  same  equation,  we  add  subscripts 
to  the  symbol  |i  to  explicitly  show  which  \x  are  the  same  and  which  are  uncorrected. 

We  use  the  symbol  t]  in  a  fashion  analogous  to  the  use  of  the  symbol  |i  to  describe  errors  unifoimly 
distributed  over  the  range  [-1,1].  We  use  these  uniformly  distributed  errors  to  account  for  worst-case 
error  limits.  As  a  general  rule,  after  we  correct  for  a  systematic  error  6,  we  then  add  errors  in  the 
simulations  that  are  uniformly  distributed  over  the  range  [-6/2,6/2]  to  account  for  the  uncertainty 
in  the  correction  we  have  applied.  We  use  smaller  ranges  in  the  simulator  only  if  we  are  able  to 
present  arguments  that  the  uncertainty  in  our  estimate  of  the  systematic  errors  is  smaller  than 
[-6/2,6/2].  We  believe  that  in  most  cases  this  a  conservative  estimate  of  the  uncertainty  in  the 
correction. 


8.1    On- wafer  TRL  Scattering  Parameter  Calibration 


These  errors  are  represented  as  transmission  matrices  that  we  cascade  onto  the  right-hand  side  of  the 
scattering  parameters  of  the  probe  head.  We  also  cascade  them  onto  the  left-hand  side  of  our 
measurements  of  the  CPW  resistors  that  we  use  to  terminate  the  signal  coming  from  the 
photodetector  (see  Fig.  6). 

8.1.1  TRL  Reference  Impedance  Offset  (from  Ref.  [4]) 

The  errors  in  determining  the  overall  magnitude  of  the  reference  impedance  of  the  TRL  calibration 
are  represented  by  the  transmission  matrix  Tp,  which  is  the  cascade  matrix  of  an  ideal  impedance 
transformer.      is  given  by 


T  = 


1 


Vi-Ar^ 


1  Ar 
Ar  1 


(11) 


where 


Ar  = 


1 

r  AC] 

AL 

[AC] 

[AC] 

2 

I  c) 

^1  ^ 

c 

+ 

I  c) 

w 

+ 

I  C) 

dc 


(12) 


The  component  of  the  systematic  error  in  the  load  method  we  use  to  calculate  the  capacitance  of  the 
CPW  is  the  worst-case  error 


AC 
C 


=  0.01 


(13) 
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The  value  in  Eq.  (13)  was  derived  from  Ref.  [5],  where  we  compared  measurements  of  CPW 
capacitance  performed  by  different  methods. 


The  component  due  to  CPW  Hne-length  variations  is 

AZ  _   1  jim 


L       20  cm 


0.000005  . 


(14) 


We  derived  this  from  a  worst-case  variation  in  Hne-length  of  ±1  jim  we  expect  on  the  masks  and 
from  the  analysis  presented  in  Ref.  [4]. 

The  component  due  to  CPW  line-width  variations  is 


C 


=  0.024  .  (15) 


w 


We  derived  this  from  field  calculations  using  the  method  of  Ref.  [6]  assuming  a  ±1  \im  error  in 
linewidths  on  the  CPW.  This  results  in  a  2  |j.m  change  in  the  gap,  which  corresponds  to  a  change  in 
capacitance  of  2.4  %.  These  calculations  are  shown  in  Appendix  1. 

We  also  add  a  component  not  discussed  in  Ref.  [4]  to  Eq.  (12)  to  account  for  errors  in  our  on-wafer 
DC  resistance  measurements.  That  factor  is 


C) 


=  V  2  — -  ,  (16) 


DC  -^DC 


where  R^^q  is  the  measured  DC  resistance  of  the  CPW  load  used  to  determine  the  capacitance  of  the 
CPW  line,  and  o^q  is  our  estimate  of  the  standard  deviation  of  the  error  in  measuring  that  DC 
resistance.  The  square  root  of  2  is  required  in  Eq.  (16)  because  the  difference  of  two  resistance 
measurements  determines  the  resistance  of  the  load,  a  measurement  of  the  DC  resistance  of  an  on- 
wafer  short  and  a  measurement  of  the  DC  resistance  of  an  on-wafer  load. 

We  found  that  our  accuracy  in  measuring  DC  resistance  depends  on  the  quality  of  the  contacts  that 
the  probe  tips  make  with  the  CPW  lines  printed  on  the  LiTaOg  wafer,  and  degrades  with  the  number 
of  connections.  For  example,  in  one  experiment,  initial  DC  resistance  measurements  were  repeatable 
to  within  about  0.01  Q,  but  after  100  contacts  on  the  same  device,  we  measured  a  standard  deviation 
of  our  on-wafer  resistance  measurements  that  increased  to  over  4  Q.  Thus,  we  estimate  the  standard 
deviation  Ojx:  of  our  DC  resistance  measurements  at  the  time  at  which  we  performed  each 
calibration.  To  do  this,  we  measure  the  DC  resistance  of  each  standard  in  both  our  initial  and  repeat 
measurements.  Then,  we  estimate  the  standard  deviation  o^c  from  the  differences  Ai?,  we  measured 
under  the  assumption  that  the  mean  of  the  A/?,  was  0.  We  use  the  formula 
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o 


DC 


2  n 


DC 


(17) 


where  h^q  is  the  number  of  measurements.  The  factor  of  2  is  required  in  (17)  because  the  A/?, 
correspond  to  the  difference  of  two  resistance  measurements,  and  we  wish  to  estimate  the  error  in 
a  single  resistance  measurement, 

8.1.2  TRL  Reference-plane  Position  Error  Due  to  Asymmetric  CPW  Short  (from  Ref.  [4]) 

We  represent  the  errors  in  the  reference-plane  position  of  the  TRL  calibration  with  the  transmission 
matrix  corresponding  to  the  scattering  matrix  5,,  which  is  an  approximation  to  the  scattering  matrix 
of  an  ideal  reference  plane  shift.     is  given  by 


where 


5.  - 


6-^ 


0  l+z6 
1+/6  0 


(1  |im)  Ti 


N  2 


(18) 


(19) 


and  where  e^,  the  approximate  relative  dielectric  constant  of  the  LiTa03  substrate  at  microwave 
frequencies,  is  43.  These  expressions  were  derived  from  approximations  in  Ref.  [4]  and  the 
assumption  of  an  error  of  ±1  |im  in  the  metal  pattern  on  the  wafer. 


8.1.3  Metal  Conductivity  (from  Ref.  [4]) 


Variations  in  metal  conductivity  and  thickness  within  a  TRL  calibration  set  also  introduce  errors  into 
the  calibration.  We  accounted  for  these  errors  with  the  simulator  described  in  Ref.  [4].  At  each 
frequency,  we  used  the  simulator  to  generate  the  scattering  parameter  matrix  5^,,  which  represents 
a  worst-case  error  in  the  on-wafer  TRL  calibration  due  to  variations  in  the  resistance  of  the  lines.  As 
inputs  to  the  program,  we  used  an  on-wafer  scattering-parameter  calibration  performed  in  the  lines 
and  introduced  a  resistance  perturbation  in  the  longest  line  based  on  DC  measurements  of  the 
resistances  of  the  line.  This  introduces  the  greatest  errors  in  the  calibration  [4]. 

Finally,  we  represent  the  errors  in  terms  of  a  transmission  matrix  corresponding  to  the  scattering 
parameter  matrix  S^'  defined  by 


0 

1 

mi2 

+ 

1 

0 

'^m22 

We  correlate  the  errors  in  Eq.  (20)  because  they  have  a  single  source,  the  resistance  of  the  line. 
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8.1.4  Metal  Conductivity  Across  the  Wafer 


Variations  in  metal  conductivity  and  thickness  across  the  wafer  also  introduce  errors  into  the 
measurements.  This  occurs  because  the  cahbration  is  performed  at  a  location  with  one  value  of  metal 
conductivity  and  thickness,  while  the  appropriate  calibration  requires  that  the  metal  conductivity  and 
thickness  be  equal  to  that  at  the  location  at  which  the  measurement  was  made.  We  account  for  these 
errors  by  performing  two  calibrations,  one  at  the  center  of  the  wafer,  where  the  lines  are  least 
resistive,  and  one  at  the  edge  of  the  wafer,  where  the  lines  are  most  resistive.  Then  we  use  the 
calibration  comparison  method  of  Ref.  [7]  to  quantify  the  differences  in  the  two  calibrations. 

We  represented  the  errors  in  terms  of  a  transmission  matrix  corresponding  to  the  scattering  parameter 
matrix  5^"  defined  by 


0 

1 

'^mll 

mlz 

+ 

1 

0 

^ml2 

Here  '  was  the  scattering  parameter  matrix  representing  the  differences  between  the  two 
calibrations  determined  by  the  calibration  comparison  method.  We  correlate  these  errors  because 
they  have  a  single  source,  the  local  conductivity  and/or  thickness  of  the  metal  film.  These  errors  were 
not  considered  in  the  analyses  we  performed  before  May  2004. 

8.2    Drift  Error  Measured  by  the  Calibration  Comparison  Method  of  Ref.  [7] 

To  determine  the  scattering  parameters  of  the  probe  head,  we  first  performed  a  coaxial  short-open- 
load-thru  (SOLT)  calibration.  Then  we  performed  two  on-wafer  TRL  calibrations.  We  calculate  the 
scattering  parameters  of  the  probe  head  from  the  first  on-wafer  calibration,  and  the  drift  in  the 
measurement  from  the  differences  between  the  first  and  the  second  on-wafer  calibration  using  the 
calibration  comparison  method  described  in  Ref.  [7].  Since  these  errors  are  generally  due  to  only  a 
few  sources  of  drift  in  the  analyzer,  they  are  highly  correlated.  Using  the  same  form  as  above,  we 
represent  the  drift  errors  in  terms  of  the  transmission  matrix  corresponding  to  the  worst-case  drift 
errors  from  the  scattering  parameter  matrix  A5",  where 


0  1 

+ 

1  0 

(22) 


In  (22),  A5  is  the  scattering  parameter  matrix  representing  the  worst-case  error  in  the  on-wafer  TRL 
calibration  due  to  test-set  drift  determined  by  the  calibration  comparison  method. 
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8.3  Error  in  the  DC  Scattering  Parameters 


We  also  determine  the  scattering  parameters  of  the  probe  head  and  the  reflection  coefficient  of  the 
CPW  load  from  on-wafer  DC  resistance  measurements.  We  add  a  resistance  of  liiOoc  to  our 
measurement  of  the  DC  resistance  of  the  probe  head  and,  since  the  resistance  of  the  on-wafer  CPW 
load  is  set  to  the  difference  of  the  measurement  of  the  resistances  of  the  load  and  a  short,  we  add  a 
resistance  of  \/2  \X20^^  to  the  measured  resistance  of  the  CPW  load. 

8.4    Coaxial  SOLT  Scattering  Parameter  Calibration 

In  our  early  work,  we  tested  photodetectors  with  2.4  mm  coaxial  connectors.  More  recently,  we  have 
focused  on  testing  high-bandwidth  photodetectors  with  1  mm  coaxial  connecters. 

We  estimate  errors  of  the  40  GHz  2.4  mm  coaxial  SOLT  calibrations  used  to  characterize  our 
photodetectors  with  2.4  mm  coaxial  connectors  with  the  "HP8510  Specifications  and  Performance 
Verification  Analysis  Software"  [8],  [9]  provided  by  the  manufacturer  of  the  vector  network  analyzer 
and  coaxial  calibration  kit  used.  We  estimate  errors  in  the  1 10  GHz  1  mm  coaxial  SOLT  calibrations 
used  to  characterize  high-bandwidth  photodetectors  with  the  "HP8510  Specifications  and 
Performance  Verification  Analysis  Software"  [10]. 

The  computer  programs  [8]  and  [10]  do  not  supply  information  on  how  the  errors  of  the  SOLT 
calibrations  are  correlated  with  frequency,  so  we  did  not  correlate  the  errors  of  the  SOLT  calibration 
with  frequency.  Thus,  our  Monte-Carlo  simulations  are  most  accurate  when  estimating  the 
uncertainties  of  quantities  determined  one  frequency  at  a  time,  such  as  the  power  spectral  density 
generated  by  the  detector  or  the  absolute  phase  of  its  response.  Because  we  do  not  calculate  these 
correlations,  the  error  analysis  is  not  suitable  for  determining  the  temporal  properties  of  the 
calibrated  pulses. 

In  this  program,  we  chose  the  parameters  most  applicable  to  our  equipment  and  procedures.  For  our 
2.4  mm  calibrations,  we  specified  the  HP8516A  test  set,  the  HP834X016  source,  the  HP85056A 
calibration  kit,  and  the  sliding  load  calibration  technique.  Although  we  used  cables  from  another 
manufacturer  not  supported  in  the  error  analysis  software,  we  specified  the  HP  85 133F  pair  of  short 
cables  because  we  felt  that  they  corresponded  most  closely  to  the  cables  we  used  in  our  experiments. 

The  program  generates  tables  of  what  we  treat  as  3a  errors  for  the  frequency  ranges  0.045  to  2  GHz, 
2  to  20  GHz,  20  to  36  GHz,  and  36  to  40  GHz.  We  generate  our  distributions  of  errors  in  our  Monte- 
Carlo  simulator  based  on  the  uncertainties  in  these  tables.  The  decision  to  treat  these  errors  in  this 
way  is  based  on  conversations  with  Doug  Rytting  and  Ken  Wong  of  Hewlett-Packard.  Doug  and  Ken 
were  involved  in  the  development  of  the  software  described  in  Refs.  [8]  and  [10]  and  with 
unpublished  work  used  to  investigate  the  proper  treatment  of  the  quantities  generated  by  Refs.  [8] 
and  [10].  The  values  drawn  from  Ref  [8]  are  summarized  in  Table  1.  The  quantities  denoted  with 
a  prime  in  the  table  correspond  to  uncertainties  in  the  characterization  of  the  reflection  coefficient 
of  the  detector,  rather  than  the  error  boxes. 
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Table  1. 

2.4  mm  connector  3-o  errors  reported  by  [8] 


<2GHz       2-20  GHz       20-36  GHz       >  36  GHz       Measurement  Condition 


„  a 
^11 

nil  'J  v  1  C\~^ 

i.OUUiXiU 

1  All  Csy  1  A~2 
1 .01  loXiU 

Ic  l~lc  l~r»  Ic  l~lc  \~c\ 

|5ii|~|522|~0,  |i2ll~|5i2|~0 

1.717x10-2 

1.7594x10-2 

2.8513x10-2 

2.8631x10-2 

|5nlH522hO,  |52,|H5i2l=l 

5.8x10"^ 

7.2x10-^ 

1.66x10-2 

1.76x10-2 

|5„|H522l  =  0,  |52jH5,2hl 

^<D2l'^ 

7.9x10"^ 

6.96x10-2 

0.127 

0.14 

l^uhfel-O,  |53jH5l2l=l 

8.492x10"^ 

8.952x10-' 

1.4017x10-2 

1.4058x10-2 

HP83440D 

1.92x10-2 

8.05x10-2 

0.1456 

0.1678 

^  From  table  "Sll  lower  worst-case  uncertainty  specifications". 

From  table  "S22  lower  worst-case  uncertainty  specifications". 

"  From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  m  in  dB  in  the  row  corresponding 
to  the  level  of  |S2,|  equal  to  0  dB  and  calculate  e^^  from  1-10"^°. 

^  From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  6  in  degrees  in  the  row 
corresponding  to  the  level  of  \S2^\  equal  to  0  dB  and  set  e^ji  equal  to  sin(6). 

^  From  table  "Sll  lower  RSS  uncertainty  specifications".  We  picked  values  in  the  rows  corresponding  to  reflection 
coefficient  |S„|  =  0.5  below  20  GHz,  |S„|  =  0.4  from  20-36  GHz,  and  |5„|  =  0.3  from  36-40  GHz.  Since  the  values 
of  e,,'  and  e^ii'  depend  on  the  reflection  coefficient  of  the  detector,  values  used  for  other  detectors  may  be  different. 

^  Values  correspond  to  preceding  row,  but  we  take  the  table  value  6  in  degrees  and  set  e^i,  equal  to  sin(0). 
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We  also  chose  parameters  applicable  to  our  equipment  for  our  1 10  GHz  1  mm  coaxial  SOLT 
calibrations.  In  the  program  [10]  for  the  frequency  range  of  45  MHz  to  50  GHz,  we  specified  the 
HPE7342  test  set,  the  HP8365xXF  source,  the  HP85059A  calibration  kit,  and  the  broadband  load 
calibration  technique.  For  the  frequency  range  of  50  GHz  to  110  GHz,  we  specified  the 
HPE7352A1 10  test  set,  the  HP836xxXF  source,  the  HP85059A1 10  calibration  kit,  and  the  offset 
short  calibration  technique.  We  also  specified  the  DIRECTXF  test  port  cables  (no  cables)  because 
no  specification  for  1  mm  test  cables  was  available  in  the  analysis. 

The  program  generated  tables  of  what  we  treated  as  3a  errors  for  the  frequency  ranges  2  to  18  GHz, 
18  to  40  GHz,  40  to  50  GHz,  50  to  75  GHz,  75  to  85  GHz,  85  to  100  GHz,  and  100  to  1 10  GHz.  We 
generated  our  distributions  of  errors  in  the  simulator  based  on  the  uncertainties  in  these  tables.  The 
values  drawn  from  [10]  are  summarized  in  Tables  2  and  3.  Here  again,  the  quantities  denoted  with 
a  prime  in  the  tables  correspond  to  uncertainties  in  the  characterization  of  the  reflection  coefficient 
of  the  detector,  rather  than  the  error  boxes.  Since  the  program  did  not  report  errors  below  2  GHz,  we 
used  the  values  reported  by  the  program  from  2  to  18  GHz  below  2  GHz. 

To  account  for  the  effects  of  bending  on  the  cables,  we  measured  the  scattering  parameters  of  a 
coaxial  "thru"  connection  after  each  on-wafer  experiment.  Since  we  know  that  this  thru  connection 
is  reciprocal,  we  averaged  the  phase  of  the  forward  and  backward  transmission  coefficients,  and  fit 
the  result  to  a  line  with  an  intercept  at  0  to  determine  the  phase  slope  C,  for  each  measurement.  We 
then  estimated  the  standard  deviation  of  the  phase  slope  under  the  assumption  that  the  mean  phase 
slope  was  0.  We  used  the  formula 


o,  = 


_!   (23) 


2n 


C 


where  was  the  number  of  measurements.  The  factor  of  2  in  the  denominator  of  Eq.  (23)  accounts 
for  the  fact  that  we  measure  the  total  phase  change  in  the  two  cables,  whereas  the  error  in  the 
measurement  is  caused  by  only  the  cables  attached  to  the  probe  head  we  are  characterizing. 
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Table  2. 

1  mm  connector  3-o  errors  reported  in  [10]  for  frequencies  less  than  50  GHz 


<2GHz       2-18  GHz       18-40  GHz      40-50  GHz      Measurement  Condition 


p  a 

0  032623 

V/ .  Vy    X*  V/ X*  ^ 

0  Q'ilGTi 

\J  t,\J ^ ^ 

0  051 144 

0  064121 

K  I^K  1=0  K  l=K  1  =  0 
Pill  P22r^'  P21I  Pi2n^ 

0.077705 

0.077705 

0.122345 

0.153653 

\SxxHS22\-O,  \s,,\-\s,,\^i 

3.235x10"' 

3.235x10-2 

5.118x10-2 

8.085x10-2 

\Sxx\-\S22\-O,  \S,MSx2\-^ 

e^2x 

3.41x10-2 

3.41x10-2 

5.955x10-2 

9.1x10-2 

\SnHS22\-O,  \s,,Hs,,\^i 

en" 

4.946x10  2 

3.9708x10-2 

6.3261x10-2 

6.8597x10-2 

NEL  detector 

5.685x10-2 

5.863x10-2 

9.887x10-2 

0.1474 

(1^22^0,  \S,,\^\S,,\^0) 

0.031688 

0.031688 

5.0441x10-2 

6.3491x10-2 

U2T  detector 
T2.D25.515.B6.116 

if 

0.3188 

0.3188 

0.2606 

0.3273 

(1522^0,  |52i|H5l2hO) 

*  From  table  "Sll  lower  worst-case  uncertainty  specifications". 

From  table  "S22  lower  worst-case  uncertainty  specifications". 

"  From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  m  in  dB  in  the  row  corresponding 
to  the  level  of  [Sjil  equal  to  0  dB  and  calculate  ^21  from  1  -  10""^°. 

^  From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  0  in  degrees  in  the  row 
corresponding  to  the  level  of  |52i|  equal  to  0  dB  and  set  e^2\  equal  to  sin(0). 

^  From  table  "Sll  lower  RSS  uncertainty  specifications"  generated  by  [10].  We  picked  values  in  the  rows 
corresponding  to  reflection  coefficient  |S„|  =  0.9  below  2  GHz,  |5„|  =  0.7  from  2  to  40  GHz,  |5„|  =  0.5  from  40  to  50 
GHz,  |S„|  =  0.4  from  50  to  85  GHz,  and       =  0.3  from  85  to  1 10  GHz.  Since  the  values  of  e,,'  and  e^,,'  depend  on 
the  reflection  coefficient  of  the  detector,  values  used  for  other  detectors  may  be  different. 

'  Values  correspond  to  preceding  row,  but  we  take  the  table  value  6  in  degrees  and  set  e^,,  equal  to  sin(6). 

« From  table  "Sll  lower  RSS  uncertainty  specifications"  generated  by  [10].  We  picked  values  in  the  rows 
corresponding  to  reflection  coefficient  |5„|  =  0.1  below  18  GHz,  |5„|  =  0.2  from  18  to  50  GHz,  |5,,|  =  0.5  from  50  to 
85  GHz,  and  |5,,|  =  0.8  from  85  to  1 10  GHz.  Since  the  values  of  e,,'  and  e^„'  depend  on  the  reflection  coefficient  of 
the  detector,  values  used  for  other  detectors  may  be  different. 
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Table  3. 

1  mm  connector  3-o  errors  reported  in  [10]  for  frequencies  greater  than  50  GHz 


i 


50-75  GHz 

75-85  GHz 

85-100  GHz 

100-1 10  GHz 

Measurement  Condition 

en 

0.042998 

0.043213 

0.053521 

0.053521 

l'5„|=|522|  =  0,  |S2,|=|5,2|  =  0 

1       ill        1       ^£,1              '1       ^11        1  1^1 

0.084111 

0.084266 

0.104882 

0.104882 

^21^ 

3.891x10"^ 

4.343x10"' 

5.421x10"' 

5.421x10"' 

4.852x10"^ 

5.924x10"' 

7.277x10"' 

7.451x10"' 

|5iilH522l=0,  \S,MSn\-^ 

e^r 

4.306x10-2 

4.308x10  ' 

5.1644x10  ' 

5.1644x10  ' 

NEL  detector 

0.1178 

0.1239 

0.1908 

0.1926 

(|522|=0,  |52.JH5.2|=0) 

en" 

4.511x10"^ 

4.514x10' 

6.6559x10"' 

6.6559x10"' 

U2T  detector 
T2.D25.515.B6.116 

0.1003 

0.1064 

0.1019 

0.1037 

(|522h0,  |52jH5,2hO) 

^  From  table  "S 1 1  lower  worst-case  uncertainty  specifications". 

From  table  "S22  lower  worst-case  uncertainty  specifications". 

'  From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  m  in  dB  in  the  row  corresponding 
to  the  level  of  |52,|  equal  to  0  dB  and  calculate  e^^  from  1-10'"^°. 

From  table  "S21  lower  RSS  uncertainty  specifications".  We  take  the  table  value  0  in  degrees  in  the  row 
corresponding  to  the  level  of  IS21I  equal  to  0  dB  and  set  e^2\  equal  to  sin(6). 

"  From  table  "Sll  lower  RSS  uncertainty  specifications"  generated  by  [10].  We  picked  values  in  the  rows 
corresponding  to  reflection  coefficient  |5,,|  =  0.9  below  2  GHz,       =  0.7  from  2  to  40  GHz,  |5,,1  =  0.5  from  40  to  50 
GHz,  |5,,|  =  0.4  from  50  to  85  GHz,  and  |5,i|  =  0.3  from  85  to  110  GHz.  Since  the  values  of  e,,'  and  e^,,'  depend  on 
the  reflection  coefficient  of  the  detector,  values  used  for  other  detectors  may  be  different. 

^  Values  correspond  to  preceding  row,  but  we  take  the  table  value  0  in  degrees  and  set  e^,,  equal  to  sin(0). 

^  From  table  "Sll  lower  RSS  uncertainty  specifications"  generated  by  [10].  We  picked  values  in  the  rows 
corresponding  to  reflection  coefficient  |5,,|  =  0.1  below  18  GHz,       =  0.2  from  18  to  50  GHz,       =  0.5  from  50  to 
85  GHz,  and       =  0.8  from  85  to  1 10  GHz.  Since  the  values  of  e,,'  and  e<(,,,'  depend  on  the  reflection  coefficient  of 
the  detector,  values  used  for  other  detectors  may  be  different. 
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8.4.1  Probe-head  Characterization 


We  incorporated  the  effect  of  the  coaxial  SOLT  caHbration  on  the  scattering  parameters  of  the  probe 
head  by  cascading  the  transmission  matrix  corresponding  to  5solt  on  the  left  side  of  the  probe  head. 
We  developed  expressions  for  the  matrix  based  on  the  relations  (26)  to  (29)  in  Ref .  [11]  between  the 
12-term  and  8-term  error  model.  We  deteniiined  5solt  from 


'SOLT 


11 


3.41 


^ 


3.41 


•(24) 


In  Eq.  (24)  the  factors  of  3  in  the  off-diagonal  terms  convert  the  3-a  limits  from  [8]  and  [10]  and 
summarized  in  Tables  1,2,  and  3  into  l-o  values.  The  square  roots  in  the  off-diagonal  terms  account 
for  the  fact  that  the  transmission  errors  from  the  SOLT  calibration  are  evenly  distributed  between 
the  error  boxes  of  the  probe  we  are  trying  to  characterize  on  port  1  of  the  network  analyzer  and  the 
error  box  of  the  second  probe  on  port  two  of  the  network  analyzer.  The  two  off-diagonal  terms  are 
equal  because  the  2-tier  algorithm  we  use  forces  the  forward  and  reverse  transmission  coefficients 
to  be  equal,  which  is  justified  based  on  reciprocity  arguments. 


We  determined  the  first  diagonal  term  in  Eq.  (24)  from  the  values  of  e,,  in  the  tables,  the 
uncertainties  in  measuring  reflection  coefficients  under  the  conditions  =  0  and  \S2\\  =  0.  The  error 
in  the  measurement  of  the  reflection  coefficient  5i,  is  equal  to  the  first  diagonal  term  in  Eq.  (24). 

Under  the  conditions  |5n|  =  |522|  =  0  and  |52i|  =  |52i|  =  1,  the  total  error  described  by  the 
measurement  of  the  reflection  coefficient  ^22  has  two  major  components.  The  first  is  the  error 
described  by  e,,  due  to  the  finite  port  match  on  port  2,  and  the  second  is  due  to  the  second  diagonal 
term  in  Eq.  (24).  Thus  we  subtracted  in  a  root-mean-square  sense  the  partial  error  e^^  we  were  not 
interested  in  from  the  total  error  ^22  to  arrive  at  the  second  diagonal  term  in  Eq.  (24)  we  were  trying 
to  determine.  We  used  the  assumption  that  the  two  error  sources  were  uncorrected  to  arrive  at  this 
expression. 


Finally,  the  factor  of  3.41  in  the  diagonal  elements  in  Eq.  (24)  convert  the  3-a  magnitude  limit  R 
from  [8]  and  the  table  into  l-o  values  for  the  real  and  imaginary  part  of  the  reflection  coefficients. 
To  arrive  at  the  factor  of  3.41,  we  assumed  Gaussian  distributions  in  x  and  y  with  variance  o.  This 
results  in  a  Gaussian  distribution  in  magnitude  r  with  no  angular  dependence.  We  then  integrated 
this  magnitude  distribution  over  a  disk  of  radius  R,  and  arrive  at  the  99.7  percentile  point  for  R  = 
3.41  a.  The  calculations  are  outlined  in  Appendix  2. 
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8.4.2  Photodetector  Reflection  Coefficient 


Because  the  detector  has  a  large  reflection  coefficient  at  its  coaxial  port,  we  treated  the  magnitude 
errors  described  by  e^/  and  the  phase  errors  described  by  e^,,'  separately.  In  the  simulations  we 
determined  the  reflection  coefficient  F^'  at  the  coaxial  port  of  the  photodetector  by  perturbing  the 
nominal  value  of  F^j  with 


In  Eq.  (25)  the  factors  of  3  convert  the  3-o  limits  from  [8]  as  summarized  in  the  tables  into  l-o 
values.  The  factor  of  Fy|F(j|  ensures  that  the  magnitude  errors  described  by  e,  i '  change  the  magnitude 
of  Fj  and  the  phase  errors  described  by  e^^^'  change  the  phase  of  F^. 

In  order  to  avoid  damaging  the  detectors  with  our  ohmmeter,  we  estimated  the  DC  reflection 
coefficient  of  the  detector  from  the  low-frequency  limit  of  its  reflection  coefficient  F^'. 

8.5    Finite  Impulse  Response  of  the  Electro-optic  Sampling  System 

The  impulse  response  of  the  electrooptic  sampling  system  is  broadened  by  several  distinct  physical 
mechanisms  [2].  Thus  the  measured  voltage  waveform  on  the  CPW  wafer  is  a  convolution  of  the 
actual  voltage  present  there  and  the  finite  impulse  response  of  the  electrooptic  sampling  system.  This 
convolution  in  the  time  domain  can  be  represented  as  a  multiplication  in  the  frequency  domain, 
which  is  how  we  choose  to  treat  these  errors.  This  section  describes  how  we  accounted  for  these 
systematic  errors,  which  we  discussed  in  greater  detail  in  Ref.  [2],  in  our  simulations. 

8.5.1  Optical  Reflection  from  the  Back  Surface  of  the  LiTaOj  Wafer 

A  small  portion  of  the  optical  sampling  beam  bounces  off  of  the  back  side  of  the  wafer,  returns  to 
the  surface  where  it  re-samples  the  electrical  signal,  and  then  re-reflects  to  the  detectors.  We  used 
the  analysis  described  in  Ref.  [2]  to  estimate  the  worst-case  magnitude  error  in  decibels  and  the 
phase  error  0^  in  degrees  of  the  measurement  system  at  each  frequency  due  to  the  first  three  optical 
round-trip  reflections  in  the  wafer. 

To  perform  the  calculations,  we  assumed  that  the  anti -reflective  coating  on  the  back  of  the  wafer 
reflected  0.05  %  of  the  incident  optical  power.  This  worst-case  estimate  was  based  on  a  reflectance 
curve  supplied  by  the  manufacturer  of  the  wafer,  and  a  measurement  we  performed  of  the  spectrum 
of  the  optical  beam,  which  showed  that  the  optical  power  was  concentrated  between  1500  nm  and 
1600  nm.  These  measurements  and  the  curve  from  the  manufacturer  are  shown  in  Appendix  3.  We 
calculated  the  reflection  coefficient  of  the  front  surface  of  the  wafer  from  the  index  of  refraction  n 
of  the  LiTa03,  for  which  we  used  a  value  of  2.12  [12]. 

We  do  not  know  if  the  optical  reflections  in  the  wafer  add  or  subtract  from  the  magnitude  and  phase 


(25) 
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of  the  final  result,  so  we  did  not  correct  for  for  this  systematic  error.  However,  we  expect  these  errors 
to  be  correlated,  so  the  simulator  adds  random  magnitude  errors  r|  to  the  magnitude  calculation, 
which  was  expressed  in  decibels,  and  random  phase  errors  r|  0^  to  the  phase  calculation,  which  was 
expressed  in  degrees  to  account  for  the  uncertainty  of  the  measurements. 

8.5.2  Radius  of  the  Optical  Beam  Traversing  the  Wafer 

We  focus  the  optical  sampling  beam  to  a  small  spot  on  the  wafer.  The  finite  radius  of  the  spot 
broadens  the  measured  pulse  slightly.  We  used  the  analysis  introduced  in  Ref.  [13],  which  assumes 
a  lowest-order  cylindrical  Gaussian  mode  of  radius  r,  to  estimate  the  error  in  the  magnitude  response 
of  the  electrooptic  sampling  system. 

We  first  estimated  that  the  radius  of  the  spot  was  about  3  [im  at  the  surface  of  the  wafer  by  passing 
the  beam  through  small  holes  in  the  metallization  on  the  wafer. 

By  Gaussian  optics,  we  estimated  the  radius  r  of  the  beam  in  the  substrate  as  r{z)  =  r(0)  (1  +  {zlz^^Y^, 
where  z  is  the  depth  in  the  substrate,  Zr=  ti  (r(0))^  niX,  n  =  2.12  is  the  index  of  refraction  of  the 
LiTaOj,  and  A,  the  optical  wavelength,  is  about  1550  nm.  Under  the  assumption  that  r(0)  is  equal 
to  3  |j,m,  the  calculation  predicts  a  beam  radius  r'  =  r(35  jim)  =  4  |j,m  at  a  point  35  |im  into  the 
substrate,  which  corresponds  roughly  to  the  50  %  point  of  the  electric  field  in  the  substrate. 

Based  on  these  estimates,  we  adopted  a  correction  for  our  magnitude  response  due  to  an  error 
calculated  using  formula  (10)  in  Ref.  [13]  with  r'  =4  \im,  our  best  estimate  of  the  average  beam 
radius  in  the  substrate.  (Note  that  the  authors  of  Ref.  [13]  use  w  to  express  the  radius  of  the  beam 
rather  than  r'.)  We  expressed  the  correction  in  decibels.  We  also  added  random  errors  VaT)  E^  to  the 
magnitude  calculation,  which  was  also  expressed  in  decibels. 

8.5.3  Finite  Temporal  Width  of  the  Optical  Pulses  Measured  by  the  Autocorrelator 

The  shape  and  temporal  width  of  the  optical  pulses  emanating  from  the  laser  do  not  affect  the  phase 
response  of  the  measurement  system,  but  do  lower  its  magnitude  response.  We  used  optical 
autocorrelator  measurements  to  correct  for  this  systematic  measurement  error. 

To  estimate  these  errors,  we  took  a  direct  Fourier  transform  of  the  autocorrelator  measurements. 
After  determining  the  error  E^  in  our  autocorrelator  measurements,  we  corrected  our  magnitude 
response  by  E^  decibels  to  account  for  this  systematic  error.  Finally,  the  simulator  added  random 
magnitude  errors  of  V2r\  E^  to  the  magnitude  calculation,  which  was  also  expressed  in  decibels,  to 
account  for  the  uncertainty  in  our  correction. 
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8.5.4  Finite  Time  the  Optical  Pulses  Spend  Traversing  the  Electric  Field  of  the  CPW  Mode 

The  optical  pulses  spend  a  finite  time  traversing  the  electric  field  of  the  coplanar  waveguide,  which 
penetrates  below  the  surface  of  the  LiTaOj  wafer.  As  discussed  in  Ref.  [13],  this  not  only  broadens 
the  impulse  response  of  the  measurement  system,  but  shifts  its  phase  response. 

To  correct  for  this  systematic  measurement  error,  we  performed  calculations  of  the  electric-field 
profile  in  the  substrate  with  the  full-wave  simulator  described  in  Ref  [14].  We  used  315  modes  in 
the  electromagnetic  calculation,  but  found  that  reducing  the  number  of  modes  to  half  that  many 
resulted  in  no  discemable  difference  in  the  calculations.  From  this  experiment,  we  concluded  that 
the  calculation  method,  at  least,  did  not  introduce  significant  errors  into  the  estimation  of  the  errors. 

The  ±2  \xm  worst-case  width  error  we  expect  in  the  40  jim  gap  of  the  CPW  will  introduce  a  ±  2 
|im/40  |im,  or  ±5  %,  error  in  our  estimate  of  the  extent  to  which  the  fields  penetrate  into  the 
substrate.  Calculations  of  the  correction  as  a  function  of  substrate  thickness  showed  that  a  ±  25  [im 
error  in  the  substrate  thickness  leads  to  a  negligible  ±0.02  %  deviation  in  magnitude  and  an  even 
smaller  ±0.002  %  deviation  in  phase. 

However,  the  anisotropic  permitivity  is  a  significant  source  of  error  that  we  did  not  attempt  to 
account  for.  The  relative  dielectric  constant  in  the  propagation  direction  (x)  and  into  the  substrate 
(y)  is  about  43,  while  the  relative  dielectric  constant  in  the  z  direction  transverse  to  the  propagation 
direction  and  parallel  to  the  substrate  surface  is  about  53  [12].  The  method  of  Ref.  [14]  cannot 
account  for  anisotropic  dielectric  constant,  but  we  expect  that  the  unequal  dielectric  constants  in  the 
y  and  z  directions  will  distort  the  fields  by  an  amount  roughly  equal  to  the  change  in  the  dielectric 
constant  in  the  two  directions,  or  about  26  %.  So  we  estimated  that  our  worst  case  error  in  the  overall 
estimate  was  ±0.26,  and  note  that  our  uncertainty  in  this  correction  could  be  reduced  significantly 
if  we  were  to  use  an  electromagnetic  field  simulator  that  could  account  correctly  for  the  anisotropic 
substrate  dielectric  constant. 

After  calculating  the  magnitude  error  E^pw  decibels  and  phase  error  Qqpw  degrees  from  the  field 
calculations  based  on  the  method  of  Ref.  [14]  and  the  formulas  in  Ref  [13],  we  corrected  both  our 
magnitude  and  phase  responses  to  account  for  this  systematic  error.  Finally,  we  added  random 
magnitude  errors  of  0.26  r|  E^p^  to  the  magnitude  calculation,  which  was  expressed  in  decibels,  and 
random  phase  errors  of  0.26  r|  O^pw  to  the  phase  calculation,  which  was  expressed  in  degrees.  Here 
we  left  these  two  errors  completely  correlated,  as  we  would  expect  them  to  be  in  practice. 


The  authors  thank  Robert  Judish  for  his  clear  and  thoughtful  explanations  of  the  statistical  principles 
and  analysis  employed  in  this  work,  and  Doug  Rytting  and  Ken  Wong  for  their  explanations  of  the 
HP  SOLT  error  analysis  software. 
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Appendix  1. 


A  summary  of  calculations  found  in  page  136  of  Dylan  Williams'  laboratory  notebook  #4.  A  change 
of  ±1  |im  in  a  metal  edge  changes  the  gap  width  from  40  \xm  to  38  |im,  yielding  the  24  %  relative 
change  in  capacitance. 
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Geometry  <ai1crons>: 
w  -  S.OOODfOl  s 
wg  -  1.650[H02 


4.000&f01 


S.OOOD-01 


Material  Parameters: 

substrate:  epsCr)  -     43.0  tan(de)  -  O.OOOo+OO 

conductivity  of  metallizations:       kappa  =    1. 0000+20 


Frequency:  f  =    l.OOOo-Ol  GHz 


Line  Elements: 
i  nductance  L 
capacitctnce  C 
resistance  R 
conductance  G 


5.291D-07  Henry/m 
4.541D-10  Farad/m 
8.408D-05  Ohm/m 
0.000t>+00  Siemens/m 


Cmho/m) 


1.06978D-08  db/irro 
1.23163D-06  1/m 


Propagation  constants: 

er(eff)  =  2.15930O+01  alpha  - 

beta  -  9.739020+00  1/m  alpha  - 

beta/beta (0)  =  4. 646820+00 

Characteristic  impedance: 

RealCZ)  »  3.413420+01  Ohm  Imag(2)  —4.316720-06  Ohm 
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Geoaetry  Calcrons): 
w    -    3.1000+01     s  -    3. 8000+01 
wg  =    1. 6600+02 


5.000D-01 


Material  Parameters: 

substrate:  eps(r)  =     43.0  tan(de)  = 

conductivity  of  metallizations:       kappa  = 


O.OOOo+OO 
1. 0000+20 


Frequency:  f  =    1. 0000-01  GHz 


Line  Elements: 
Inductance    L  > 
capacitance  c  = 
resistance   R  = 
conductance  G  = 


5.1630-07  Henry/m 
4.6530-10  Farad/m 
8.2860-05  Ohm/m 
O.OOOO+OO  Siemens/m 


Cmho/m) 


Propagation  Constants: 
erCeff)  =  2.15893O+01 

beta  =  9. 738200+00  1/m 

beta/betaCO)  -  4.64643O+00 


al pha  = 
alpha 


1.080270-08  db/mm 
1.24371D-06  1/m 


Characteristic  impedance: 

ReaKz)         -  3. 331200+01  Otm     InagCZ)  —4. 254420-06  ohm 
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Appendix  2. 


We  develop  the  factor  of  3.41  in  the  diagonal  elements  in  Eq.  (24)  that  convert  the  3-a  magnitude 
limit  R  from  [8]  and  the  table  into  l-o  values  for  the  real  and  imaginary  part  of  the  reflection 
coefficients.  To  arrive  at  the  factor  of  3.41,  we  assumed  Gaussian  distributions  in  x  and  y  with 
variance  a=l  with  functional  form 


271  271 


(26) 


where  /  is  the  distribution  of  errors,  x  and  >-  correspond  to  the  real  and  imaginary  coordinates,  r 
corresponds  to  the  distance  from  the  origin.  If  F(i?)  if  the  probability  that/is  inside  a  disk  of  radius 
R, 


R  2n 


0  0 


(27) 


where  R  =  sJ-2\ii{\-z).  We  thus  arrive  at  the  following  table  relating  1-,  2-,  and  3-o  limits  and  R. 


Table  2.1  Relations  between  1-a.  2-a.  and  3-o  limits  and  R. 


Limit 

1-  o 

2-  a 

3-  a 


Z(%) 
68 
95 

99.7 


R 
1.51 
2.45 

3.41 


From  the  table  we  arrive  at  the  99.7  percentile  point  fori?  =  3.41a.  This  is  close  to  the  value  of  three 
we  might  have  anticipated,  showing  that  the  x-y  distribution  differs  only  slightly  from  the  circular 
distribution,  and  the  difference  in  distributions  is  not  very  important. 
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Appendix  3. 


Measurements  of  the  optical  spectrum  (see  page  102  of  Dylan  Williams'  laboratory  notebook  #4) 
and  a  reflectance  curve  for  the  optical  antireflecti  ve  coating  on  the  back  of  the  LiTaOj  wafer  supplied 
by  the  manufacturer. 
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Data  recieved  from  the  manufacturer  characterizing  the  antireflective  optical  coating  they  deposited 
on  the  back  of  the  LiTaOj  wafers. 
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